Quantum theory of spontaneous and stimulated emission of surface plasmons 
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We introduce a quantization scheme that can be applied to surface waves propagating along 
a plane interface. An important result is the derivation of the energy of the surface wave for 
dispersive non- lossy media without invoking any specific model for the dielectric constant. Working 
in Coulomb's gauge, we use a modal representation of the fields. Each mode can be associated with 
a quantum harmonic oscillator. We have applied the formalism to derive quantum-mechanically 
the spontaneous emission rate of surface plasmon by a two-level system. The result is in very good 
agreement with Green's tensor approach in the non-lossy case. Green's approach allows also to 
account for losses, so that the limitations of a quantum approach of surface plasmons are clearly 
defined. Finally, the issue of stimulated versus spontaneous emission has been addressed. Because 
of the increasing density of states near the asymptote of the dispersion relation, it is quantitatively 
shown that the stimulated emission probability is too small to obtain gain in this frequency region. 



I. INTRODUCTION 

Quantum theory of light is a useful tool to describe 
microscopic interactions between light and matter. The 
electromagnetic state is represented by photon num- 
ber states and the electromagnetic field becomes an 
operatoi^. Such a description of light provides a quanti- 
tative description of absorption, spontaneous and stim- 
ulated emission of photons by a two-level system. In 
particular, it allows to derive a quantitative treatment 
of light amplification. It also predicts pure quantum ef- 
fects, such as photon coalescence or antibunching. Quan- 
tum theory of light can be extended to non-dispersive 
and non-lossy media. Each photon in the material cor- 
responds to the excitation of a mode characterized by 
a wave vector k and circular frequency cj, such as /c = 
nuj I c, where n is the refractive index of the medium and 
c the light velocity in a vacuum. It is the purpose of this 
paper to introduce a quantification scheme for surface 
waves propagating along an interface. 

It is well known that electromagnetic surface waves 
called surface plasmons exist at interfaces between metals 
and dielectrics . Their quantum nature has been demon- 
strated by energy loss spectroscopy experiments on thin 
metallic films reported by Powell and SwaiP. Single op- 
tical plasmons have been excited recently along a metal- 
lic nanowire^ ^ Surface plasmons are associated with col- 
lective oscillation of free electrons in the metal at the 
surface. Similar electromagnetic fields exist also on po- 
lar materials and are called surface phonon-polariton. 
Both surface plasmon-polaritons and surface phonon- 
polaritons propagate along the interface and decrease in 
the direction perpendicular to the surface. Such a reso- 
nance is therefore called surface wave in a more general 
way. Most studies deal with a plane interface between air 
or vacuum and a non-lossy material. In this case, it is 
well knownFthat a surface wave can exist if the dielectric 



constant e(uS) has a real part lower than —1. 

Losses are often a serious limitation for many practi- 
cal applications envisionned for surface plasmons. This 
problem could be circumvented by introducing gain in 
the system. Studies have been made in such a way with 
metallic nanoparticles embedded in a gain medium both 
numerically with dye molecules^ or quantum dots^ ^ and 
experimentall}^. Seidel et al. reported the first exper- 
iment demonstrating the amplification of surface plas- 
mons on a flat silver film surrounded by a solution of dye 
molecules^^. Since then, a few studies have dealt with 
stimulated emission of s urface plasmons on fiat inter- 
faces both experimentall}^^^ and theoreticall \P^ . Such 
works have pay ed the way to active plasmonica^^H^^ and 
nanolasers^^l^ or more precisely to spaser^l^, or sur- 
face plasmon amplification by stimulated emission intro- 
duced by Bergman and Stockman and demonstrated ex- 
perimentally recentl}/^^ 

It is clear that a quantum treatment of surface plas- 
mon could be useful for many applications. For instance 
an efficient single photon emitter could be optimized^. 
A quantum treatment allows to model stimulated emis- 
sion and therefore to specify gain conditions and laser 
operation. It could also allow to analyse pure quan- 
tum effects for surface plasmons such as single plasmon 
interferences, quantum correlation^^, bunching, strong 
coupling regime^^ or single photon excitation of sur- 
face plasmon^^ — . To our knowledge, the first quanti- 
zation scheme for surface plasmon on a metallic surface 
has been reported by Elson and Ritchi^^. In their work, 
the metal is characterized by a non-lossy Drude model 
so that real optical properties cannot be included. Using 
Green's approach, Gruner and Welsch introduce a quan- 
tization scheme for electromagnetic fields in dispersive 
and absorptive materials^^. It should hence be possi- 
ble to quantize the field associated with surface waves 
using their model. Note that due to losses, they can- 
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not obtain operators for modes but only local operators 
: one recovers the usual creation/annihilation operators 
in the limit of zero losses. A related work, reported in 
the early nineties by Babiker et al., dealt with the quan- 
tization of interface optical phonons in quantum well, 
which could appear also as a confined surface phonon in 
a heterostructure^^ . 

In this paper, a quantization scheme that is not based 
on a specific model of the dielectric constant is intro- 
duced. The aim is to quantize the field by accounting for 
the experimental dispersion properties of the medium. 
The procedure follows the quantization scheme for pho- 
tons in a vacuum. We will first introduce a classical mode 
description of the surface waves and discuss the disper- 
sion relation. A key issue for quantization is the defini- 
tion of the energy of surface waves for dispersive lossy 
media. The problem of electromagnetic energy in a dis- 
persive and lossy medium has been recently addressed in 
a paper by Stallinga^^. The third section addresses the 
problem of the electromagnetic energy associated with 
surface waves in a simpler case following Landau and 
Lifchitz for non- lossy dielectric material^^. The quanti- 
zation scheme is finally described in the fourth section. 
In order to check our results, we apply our formalism 
in the fifth section to the calculation of the spontaneous 
emission of a two-level system in the presence of surface 
plasmons. The Purcell factor (i.e. the local density of 
states normalized by the vacuum density of states) and 
Einstein's coefficients are also derived using this model. 



II. MODAL DESCRIPTION OF SURFACE 
WAVES 
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FIG. 1: Surface wave on a plane interface. The surface mode 
is characterized by its circular frequency uj and the projec- 
tion of the wave vector along the interface K. K and z are 
unit vectors along and perpendicular to the plane interface 
respectively. 



Let us consider surface waves propagating on a plane 
interface at 2: = separating two semi-infinite media 



(Fig. [T]). One of them is a vacuum or air and the sec- 
ond is a metal or a polar material. A surface mode is 
characterized by its circular frequency uj and the projec- 
tion of the wave vector K on the plane perpendicular to 
the z-axis. The material has a dielectric constant e(ijj). 
We use Coulomb's gauge (divA(r, t) = 0) to write the 
magnetic and electric fields: 



B(r,t) 
E(r,t) 



VxA(r,t) 
dA{Y,t) 
dt 



(1) 
(2) 



The field produced by any distribution of sources in the 
presence of an interface can be computed using Green's 
tensor. By extracting the pole contribution, it is possible 
to derive the general form of the surface plasmon field. 
The details of this procedure can be found in Ref. |35l 
The corresponding vector potential can be cast in the 
form: 



A(r,t) 



(27r) 



2 <^kUk(^) exp(iK.r) exp(— za;5pt)+c.c. 

(3) 

where c.c. stands for complex conjugate. In this equation, 
K is a real wave vector parallel to the interface and the 
circular frequency uOgp is a complex root of the equation: 



K 




(4) 



The term q^k is an amplitude associated with wave vector 
K in the decomposition. The vectors uk(^) are given by: 



1 f ^ K 



(5) 



where L{uJsp) has the d imension of a length and will be 
fixed later by Eq. (B20) to normalize the energy of each 
mode. 7j is the projection of the wave vector along the 
z-axis, j = 1 in the region z > 0, and j = 2 in the region 
z < 0, so that = ej{uJsp)oj1p/ — - The sign of 7^ 
is then chosen such as the field goes to zero when z goes 
to ±00. Let us note that in the non-lossy case, 71 and 
72 are purely imaginary, so that the electric field decays 
exponentially along the z-axis. K and z are unit vectors 
directed along K and the z-axis respectively. 

Figure [2] shows the dispersion relation of surface plas- 
mons as well as the variation of the imaginary part of 
the frequency with on a plane interface of silver. To 
perform the calculations when a complex frequency is 
needed, it has been useful to fit the experimental values 
of the dielectric constant e(a;) given by Ref. 36 with an 
analytical model. The real part of the silver dielectric 
constant is very well represented by a Drude model given 
in Ref. 15, For the imaginary part we add to this Drude 
model a conductivity term, so that the modelized dielec- 
tric constant is in very good agreement with the experi- 
mental data. The model we used is given in appendix [A] 
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FIG. 2: Dispersion relation of a surface plasmon on a plane 
interface between air and silver (solid line, left axis) and vari- 
ation of the imaginary part of cj (dashed line, right axis). The 
dispersion relation has been obtained using the silver dielec- 
tric constant given in Appendix [A] 



In this example, lm{ujsp) is small, less than 5%, compar- 
ing to Ke{uJsp). In other words, the lifetime of the surface 
mode is long enough to have a few tens of periods for the 
oscillating electromagnetic field on the asymptotic part 
of the dispersion relation and hundreds of periods on the 
linear part, close to the light cone. Note that this point 
seems to be rather general. Indeed we found similar ra- 
tios for many other materials supporting surface waves. 

From the dispersion relation, it is possible to derive 
the density of states. To this aim, it is convenient to 
introduce a virtual box, which is in fact a virtual square 
in the x—y plane of sides Lx and Ly and size S = LxXLy. 
Born- Von Karman's conditions yields a quantized value 
of the wave vector Kx = Ux^tt/Lx^ Ky = ny2i\ jLy^ where 
'^Xj'y are relative integers. Let us note that a different 
expansion of the surface plasmon field can be used with 
a complex wave vector and a real frequency as discussed 
in Ref. [35] We stress here that the Born- von Karman 
procedure imposes a real wave vector. It follows that the 
relevant dispersion relation has no backbending as seen 
in Fig.|2j The reader is refered to Ref.[35lfor more details. 
Substituting the discrete sum ^ over the quantized 

wave vector K and the discrete amplitude SAj^ to / -^^^ 
and a-K respectively, the vector potential can be cast as 



A(r, t) = ^ ^kUk(^) exp(iK.r) exp(- 

K 



-iujt) 



c.c. 



(6) 



where we have omitted the subscript sp for the circular 
frequency uj. We can insert this form in Eqs. ([T]) and ^ 
to obtain the electric and magnetic fields. Introducing 
the notations kj = K + 7jZ and bK(^) = kj x uk(^), we 



E(r,t) = i ^y^^ujAj^uj^{z) exp(zK.r) exp(— zcjt) + ({^) 

K 

B(r,t) = i A-Khi^{z) exp(zK.r) exp(— zcjt) + c.^7h) 



K 



III. ENERGY OF A SURFACE WAVE 

The quantization procedure is based on the fact that 
the energy of the field has the structure of a sum of har- 
monic oscillators. It is thus a key issue to derive the 
energy of the surface plasmon field. In this section, we 
give a brief outline of the derivation and leave the details 
to appendix [BI In a vacuum, the energy density is given 



(8) 



The electromagnetic energy in a lossy dispersive ma- 
terial is a more subtle issue. This problem has been ad- 
dressed for the first time by Brillouin^^. He considered 
a very simple case, with two perfectly monochromatic 
waves in the material. Landau and Lifchitz analysecP^ 
the energy of an electromagnetic field in a non-lossy dis- 
persive medium, whose frequencies form a narrow con- 
tinuum around the mean frequency ujq. They dealt with 
fields such as E = Eo(t) exp(— zcjq^), Eo(t) varying slowly 
over the period 27r/uJo- In the appendix, we follow this 
method. The main idea is to derive the work done 
by an external operator to build adiabatically the field 
amplitude. This work is equal to the total amount of 
electromagnetic energy of the surface waves for a non- 
lossy medium. Note that more recently, Stallinga de- 
rived an expression of the energy for dispersive and lossy 
materials^^. The result is the same provided that e is 
replaced by Re(e). This suggests that it is possible to 
neglect losses in the calculation of the energy. Actually, 
it is essential to deal with a non-lossy medium to have 
well-defined modes. A key issue regarding this approxi- 
mation is whether the dispersion relation is modified by 
the presence of losses. Indeed, the density of states crit- 
ically depends on the dispersion relation. We compared 
the dispersion relation obtained using Re{uJsp) for a lossy 
medium with the dispersion relation with a non-lossy 
medium in the case of silver. We found a relative dif- 
ference between the two dispersion relations always less 
than 1.5 x IQ-^. 

We will thus neglect the losses of the medium in the 
derivation of the energy. The calculation outlined in ap- 
pendix E] gives the total energy of the surface waves: 



(9) 



K 



We emphasize that this convenient expression for the en- 
ergy is obtained using the right normalization condition 
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on L{ uj) o r equi valen tly on uk(^) given respectively by 
Eqs. p20l) and (|B22|). 



EMISSION RATES: COMPARISON WITH 
THE CLASSICAL CASE, EINSTEIN'S 
COEFFICIENTS 



IV. QUANTIZATION OF SURFACE WAVES 

We now turn to the quantization of the electromag- 
netic field of surface plasmons. We first notice that the 
expression e^uo'^S [Ak^k + ^k^k] of the energy for each 
mode K, has the structure of the energy of a harmonic 
oscillator, hence the quantized hamiltonian: 



H 



E 

with the equivalence 



K 



fiuo 
~2 



(10) 



2eoCo'5' 



2eoCo'5' 



(11) 

(12) 



The surface wave field is thus quantized by associa- 
tion of a quantum-mechanical harmonic oscillator to each 
mode K. We introduce aj^ and ax which are respectively 
the creation and annihilation operators for the mode K. 
As in the harmonic oscillator theory, aj^ and ax act on 
surface wave number states |nK) which are eigenvectors 
associated with eigenvalues (nx + l/2)fia; of the Hamil- 
tonian (nK is an integer). Operators (respectively 
ax) allow to create (respectively destroy) a quantum of 
energy huj according to the operating rules^: 



«kI^k) 
«k|^k) 




(13) 
(14) 



Different surface modes are independent so that their 
associated operators commute: 



(15) 



We can now write the fields as operators acting on the 
surface plasmon number quantum states \nY^\ 



K 



2^ 



uk(^) ttK exp(iK.r) exp(— zcjt) 



B(r, = ]/ 260008 ^^^^"^ exp(zK. 
where h.c. denotes the hermitian conjugate. 



.r) exp(— zcjt) 
h.c. (17) 



A. Spontaneous emission of a dipole above a 
metallic interface 



The quantization scheme that we have introduced al- 
lows to derive an expression of the electromagnetic field 
using operators. Hence, we can write interaction hamil- 
tonians and describe the coupling between light and mat- 
ter. In order to test this quantization procedure, we per- 
formed the calculation of the lifetime of a two-level sys- 
tem placed in the vicinity of a metal-vacuum interface 
so that surface plasmons can be excited. This result is 
interesting as the lifetime can also be computed using 
a classical approach as shown for instance by Ford and 
Webei^^. More specifically, they showed how to find the 
surface plasmon contribution to the lifetime by extract- 
ing the pole contribution. By comparing both results, 
we can assess the validity of the quantum theory of sur- 
face plasmon within the approximation of a dispersive 
but non-lossy medium. 



1. Quantum calculation 

In the quantum approach, we first derive the decay 
rate associated to the spontaneous emission of surface 
plasmons of a two-level quantum system close to an in- 
terface, using Fermi's golden rule. This gives the surface 
plasmon spontaneous emission rate as a function of the 
matrix element (2|D|1) = D12 of the dipole moment op- 
erator D. The details of the calculation are given in 
appendix [C] 

We obtain the following expression for the spontaneous 
emission rate: 



7spont(Dl2,^0,^) 



12 



2eoh 



K 



AK 



duj Le//(^,di2,a;o) 

(18) 



in which di2 = D12/ID12I is the (possibly complex) po- 
larization of the dipole, di2,z = di2.z, di2,// = di2 — 
di2,z^' We introduced the effective length of the surface 
plasmon mode Le//(2;, di2, cjq), 



1 



exp(2z7iz) 



^h.c. (16) Le//(z,di2,a;o) L{ujq) 



di 



2,// 1 



■ e(^o)|c^i2,2p 



It will be seen later that this length allows to define an 
effective volume of the plasmon mode. 

For comparison with the classical calculation, we nor- 
malize 7spont(Di2, cjo, ^) with the spontaneous emission 
rate of the same two-state quantum system in a vacuum, 

given h^-f%^^^ = This gives the Purcell factor 



(19) 
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associated to the emission of surface plasmons: 



-K 



2ujI duj Leff{z, di2,OJo) 



(20) 



which does not depend anymore on the amphtude of D12, 
but only on its polarization di2, its frequency ujq and 
its distance to the interface z. As expected, the Pur- 
cell factor decreases exponentially as the dipole goes far- 
ther from the interface, and can have rather high values 
(see Fig. [3] and comments below) as ujq gets closer to the 
asymptotic frequency of surface plasmons if the dipole is 
not too far from the interface. 

This Purcell factor can also be cast under the form: 



Fp(di2,cJo,^) = ^'(^o) 



Veff{z,di2,UJo) StT 



(21) 



where the (global) density of states of surface plasmons 
g{uj) is given by g{uj) = and Veff{z,di2,uJo) = 

SLeff{z,di2,uJo) is the volume of the surface plasmon 
modes of frequency uq for a dipole polarization di2 in 
which the emission occurs. Eq. (21 ) is thus similar to the 



Purcell factor Fp of a dipole interacting with a single 
damped mode^^. {Fp = Q^^^^ or equivalently Fp = 
ujg{uj)^-^ using the density of states of the single mode 

at resonance, g{uj) = f^-) 

When dealing with an isotropic distribution of dipoles, 
the average of the rate of spontaneous emission ( 18 ), over 



the orientations of the dipole D12, should be considered. 
Let us first introduce the total effective length of the sur- 
face plasmon mode, defined as the inverse of the average 



of 



Leffiz,di2,UJo) 



over the directions of di2: 



Leff, totalis J ^0) 



1 exp(2z7iz) 
3 L{ujo) 



l + |e(a;o)|], (22) 



Calculating the averaged rate of spontaneous emis- 
sion then amounts to replacing Le//(2;, di2, cjq) by 
Lef totalis, ujo) in Eq. 



normalized by its value in a vacuum. Both expressions 
give the normalized local density of states, which is a 
classical quantity. They are therefore equal, that is the 
normalized radiated power gives the normalized sponta- 
neous emission rate. The normalized radiated power can 
be expressed as a function of Green's tensor G(r, r',co') 
of the system: 



^p,d(di2,^o,^) Im 



dl2'G{zz,zz,uJo)di2 . (24) 



Following the steps detailed in Ref. ^35^, the pole con- 
tribution Ggp of Green's tensor of a plane interface can 



be derived, and inserted in Eq. (24). We use here the 



pole contribution of the surface plasmon with a complex 
frequency (see Ref. 35). The details of the calculation 
are given in appendix [D] One finds for the normalized 
radiated power in the non- lossy case: 



Fp^ci{di2,uJo,z) = -^K'^--R{K,ujo) 
UJK duo 



exp(2z7iz) 



di2,//P - e(^o)|c^i2,^P 



. (25) 



Using Eq. ( |D9[ ) and comparing Eq. ([25| to Eq. ( |2Q[ ), we 

see easily that Fp (di2 ^ojq^ z) = Fp^ci {^12 , ^0 7 ^) • We thus 
recover the quantum spontaneous emission rate in the 
non lossy limit of the above classical approach. This re- 
sult is not surprising. Indeed the normalized spontaneous 
emission rate yields the local density of states. The latter 
is a classical quantity. In the quantum approach, it has 
been calculated using the dispersion relation. In the clas- 
sical approach, it has been calculated using the Green's 
tensor. We have thus checked that the mode approach 
and the Green's formalism approach are equivalent. We 
now go one step further and compare the quantum ap- 
proach (without losses) with the Green's tensor approach 
that accounts for losses. We compute the spontaneous 
emission rate for both cases in order to assess the role of 
losses. 



"yspont,total{\^12\,^0, z) 



uJo\Di2\\ ^dK 
2eQh duo 



(23) 



Leff,total{z,UJo) 

More details are given in appendix [C) 

2. Classical approach 

In the previous section, we considered a two-level quan- 
tum system having a given polarization di2 and Bohr 
circular frequency ujq , and we normalized its spontaneous 
emission rate by its value in a vacuum. The power ra- 
diated by a classical harmonic dipole having the same 
polarization di2 and a circular frequency ujq can also be 



3. Comparison with the lossy case 

The lossy and non lossy emission rates are compared 
using Eqs. ([20| and (D8). The result is seen on Fig. [s] 
for a dipole located at three different distances of a silver 
surface (10, 75 and 250 nm). It appears that the dif- 
ferences between both curves in this case are still very 
small as long as the frequency is not too close from the 
asymptote of the dispersion relation. This is not sur- 
prising considering the fact that at this asymptotic value 
the losses are the most important. Moreover, when the 
distance between the dipole and the interface increases, 
the part of the electromagnetic field due to the higher 
surface plasmons wave vector decreases, so that the part 
due to the surface plasmons lying on the linear part of the 
dispersion relation is more important. As seen in Figj2j 
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these surface waves have less losses and the quantum ap- 
proach is thus more accurate. It follows the important 
conclusion that the non-lossy medium approximation in 
the quantum treatment is reasonable to deal with surface 
waves provided that the frequency is not too close to the 
asymptotic value. 

Note that Re f.'40^ provides an expression for the surface 
plasmon emission rate of a dipole close to an interface, 
which can be compared to ours (details not given here). 
In the non lossy case, it can be proved analytically that 
their results lead to the same normalized emission rate 



as Eq. (20). In the lossy case, one finds a normalized 
emission rate close to the one used here (Eq. (|D8|)), al- 



though they are not rigorously equal. Our method gives 
an expression of the normalized surface plasmon emission 
rate as a sum over the modes K (see Eq. (D4)), which 
provides a better understanding of the difference between 
the lossy and the non lossy cases. 
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FIG. 3: Normalized emission rate Fp of a vertical dipole lo- 
cated at 10 nm (top), 75 nm (center) and 250 nm (bottom) 
from the surface when taking into account losses (green con- 
tinuous curves) comparing to the non-lossy cases (blue dash- 
dotted curves). 



B. Einstein's coefficients 

A quantum approach for surface waves allows us also 
to derive easily Einstein's coefficients for spontaneous 
and stimulated emission. The same example of a dipole 
above the interface is taken. Once again, it is pos- 
sible to follow the approach described for photons in 
Ref. [T]for instance. Einstein's coefficient for surface plas- 
mon spontaneous emission has already been calculated: 
^21 = 7spont(Di2,^o,^) (see Eq. ([l8|). In order to ob- 
tain Einstein's coefficie nt fo r stimulated emission, one 
needs to start from Eq. ( |C3| ). In this equation, the term 
proportional to is the matrix element for stimulated 
emission. We note {W{uj)) the energy density of the ra- 
diation per unit surface and we assume that it varies 



slowly for frequencies near luq. The total energy in the 
single mode n-K is now replaced by: 



nj^huj S j duj{W{uj)) 



(26) 



The transition rate due to stimulated emission can thus 
be written: 



7stim(Dl2,^0, z) 

27r 



1 

■2"^ 



2 , dc^(ipy(a;));^|Di2.UK(r)r - uj^) (27) 



It follows that Einstein's coefficient for stimulated 
emission in mode K, B21 = 7st2m(Di2, ^0, ^)/(^(^o)), 
is given by: 



21 



7r|Di2p exp(2i7iz) 



L{ujo) 



^12,// 



cos ( 



K 
— ( 

7i 



^12,; 



(28) 



where <ii2,// and di2^z are defined above, and is the an- 
gle between the projection of D12 on the interface and K. 
When dealing with an isotropic distribution of dipoles, 
B12 should be averag ed over the directions of D12, in the 
same way as in Sec. (VAl). We get the total Einstein 



coefficient for stimulated emission in mode K 



B' 



21,total 



7r|Di2p exp(2i7iz) 



[1 - ^(^0)] 



(29) 



To describe the amplification of a surface 
plasmon beam by an amplifying medium, it 
is interesting to derive the ratio r{uJo^z) = 



^2*1(10121,^0, z)/B2i dotal 



ID 



12 



(0 



where 



A^Y(|Di2|,cJo,^) = 7spon^(^pd(di2.^o,^)) Stands 
for the total spontaneous emission rate of the dipole 
close to the interface ((z) denotes interface), and Jspont 
is given above. It can be computed with Eq. (24), using 
Green's tensor of a plane interface (this rate includes 
all the waves that can be emitted, not only surface 
plasmons). (•) stands for average over the orientations 
di2 of the dipole. r(a;o, z) gives the threshold energy per 
unit surface Wc{ojo) at which the stimulated emission 
rate equals the spontaneous one. It can be written as: 



(Fg,(di2,C^0,^)) 

exp(2z7i2;) [1 - e{uJo)] 



L{coo) (30) 



where r^{uJo) = is the ratio of the Einstein's coeffi- 
cients in a vacuum. Fig. |4] shows r(a;, z) as a function of 
u. From to approximately 35 nm, the ratio decreases, 
mainly because the (total) spontaneous emission rate 
decreases. Above 35 nm, the ratio increases, as ^21 de- 
creases, because of the exponential decay of the surface 
plasmon field away from this interface. Note that for a 
given frequency cj, r{uJo^ z) rises at lower distances for uj 
closer to the surface plasmon asymptote frequency. 
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FIG. 4: Ratio of Einstein's coefficients r{(jj, z) for (total) spon- 
taneous emission and stimulated emission of surface plasmons 
as a function of the distance to the interface, for several values 
of uj (see legend). 



These results can be used to calculate the amount of 
power that undergoes stimulated emission of surface plas- 
mons in a gain medium in close vicinity of the silver in- 
terface. We consider a parallel beam of surface plasmons. 
and suppose that they are excited via a grating or a prism 
by a He-Ne laser whose emission has a linewidth of about 
Acj = 10 MHz centered around ujq = 2 eV, and that they 
carry P = 1 mW of power per jiia. The spectral power 
at maximum, assuming a lorentzian profile, is given by 
Puj{^q) = P/ttAcj, and the associated spectral energy 
per unit surface is W{ujo) = P^{ujQ)/vg = _ ^. where 



9 ~ TvAojVg 

Vg = ^ is the group velocity of surface plasmons, close 
to c below the asymptote frequency. The spectral energy 
per unit surface at maximum of these surface plasmons 
is then W{ujo) ^ j^a^c ^ 10^ ^.nm^. This value is far 
above those of Fig. [4| stimulated emission in a freely 
propagating surface plasmon beam is several orders of 
magnitude higher than spontaneous emission. 



VI. CONCLUSION 

In this paper, we have extended previous work on 
quantization of surface plasmons by introducing a for- 
malism that can use experimental values of the dielectric 
constant instead of using a specific model for the free elec- 
tron gas. The key step is the derivation of the energy of 
a surface plasmon in a dispersive non-lossy medium. The 
standard quantization scheme in Coulomb's gauge yields 
the quantum form of the field. This scheme can be ex- 
tended in a straigthforward way to thin metallic films. To 
illustrate the formalism, we have derived the spontaneous 
emission rate of surface plasmons by a two-level system 
placed close to an interface supporting surface waves as 
well as Einstein's coefficients. This quantized theory of 
surface plasmon will be useful to analyse specific quan- 



tum effects such as antibunching, single plasmon inter- 
ference, quantum coherence properties, but also to derive 
the interactions of surface waves with other quantum ob- 
jects, as quantum wells for example. 



Appendix A: Dielectric constant model 

When cases with losses are considered in this paper, we 
must sometimes consider the case of complex frequencies. 
An analytical model for the dielectric constant is needed 
to evaluate e{uj) when u is complex. The dielectric con- 
stant model which has been used in this paper is the 
following: 



e{uj) 



(Al) 



where we take the values e^o = 5, Tiujp = 9.1 eV and 
np = 0.021 eV of Drude's model given by Ref. Ol We 
add also a conductivity term to have a better fit of 
the imaginary part of the dielectric constant, we take 
ha/ Co = 1.8 eV. Fig. [s] compares the experimental datsP^ 
and both Drude's model from Ref..l5iand the model used 
in this paper in Eq.(|Al|). 



-10 
-20 




2.0 

hu (eV) 



FIG. 5: Real and imaginary part of the dielectric constant for 
silver. The experimental data from Ref. [36] in green dotted 
line, Drude's model from Ref. TF in red and the fit used in 
this paper in blue. 



Appendix B: Derivation of the energy 

In this section, we focus on the derivation of the elec- 
tromagnetic energy associated with surface waves. The 
main idea is to derive the work done by an external 
operator to build adiabatically the field amplitude in a 
non-lossy medium. The energy balance from time t = 



8 



to t = T reads: 



U 



E 



H 



dB 

~dt 



d^r (-j.E) 

in which dt J d^r (—j.E) is the operator's work on the 
system between t = and t = T. T must be large for 
this work to be adiabatic. U then does not depend on 
T. In other words, this is also the total electromagnetic 
energy of the system. Due to the exponential decrease 
along the z-axis, the contribution of the Poynting ve ctor , 
which should appear in the left term of equation (Bl), 



drops to zero. We want hence to derive the first term of 



equation (Bl) to obtain the electromagnetic energy. 



The first step is to introduce time-dependent ampli- 
tudes in Eq. ([g]): Ak is then replaced by A-K{t) in which 
^k(^) = X f{t)' This accounts for operator's work. 
For the sake of convenience, we take a 2T-periodic func- 
tion for which /(O) = /(2T) = and /(T) = 1, so that 
we can write f{t) = fn exp (^i27rn^). T is the typi- 
cal time of variation of the amplitude Ak(^). T has to be 
taken sufficiently large to consider the work done by the 
operator as adiabatic. (For instance f{t) = sin(7r2|^).) 
Let us first derive the energy in the material medium 
denoted medium 2 for 2: < 0. We have: 



E(r,t) = -dtA{r,t) = - ^ exp(zK.r)uK jz) [-icuAK. jt) dtAj^{t)] exp{-iujt) ^ c.c. 



(B2) 



K 



Using A^iit) = Ayifit) and the Fourier series expansion of f{t) given above, the field can be cast as a sum of terms 
varying as exp[— i(cj + 



c.c. 



(B3) 



C.C. 



(B4) 



E(r,t) = ^exp(iK.r)uK(2;)AK^i {uj - /n exp i - uj^ t 

K n 

Hence the displacement vector 

D(r, t) = eo ^ exp(zK.r)uK(2;)AK ^ - Y") ~ Y") ^ ( Y" ~ ^) ^ 

K n 

and its time derivative 

^ x^/ \ / \ / \ ^ / 7rn\2 „ / 7rn\ r /7rn \ " 

dtD{Y,t) = eo2^exp(zK.r)uK(^)AK2^ V~ ^) ~ 7^7 V \ Y ~^)^ 

K n 

As 1/T <C c j, we have {uo — ej{(jj — ^ uj'^ej{uj) — ^^ du^'''^ ' Taking the inverse of the Fourier series expansions 
in Eq. (B5), we get 



c.c. 



(B5) 



dtT>{r,t) = eo^exp(zK.r)uK(^) i^u'^ej{u)Aj^{t) ^ i^^^^^-^^^ (B6) 



Hence 



E{r,t).dMv,t) = 

-eo ^ exp[i(K+K').r]uK{z).UK'{z)[-iojAK{t)+dtAK{t)]\oj'hj{oj')AK'{t) 



K,K' 



.11. (, .l\ 



-eo exp[z(K-K0.r]uK(^).u;,,(^)[-ic^AK(t)+5,AK(t)] [j\,{yj')A\^\i)-i^^^^^ 

+ C.C. (B7) 

Integrating this term over the surface S and using J dx J d?/exp(i(K — K').r) = S'(5k,k' is the Kroenecker symbol. 
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which verifies ^k,k' = 1 if K = K', ^k,k' = else), we find 



dx j d?/E(r,t).5^D(r,t) 



c.c. (B8) 



The first term on the right hand side of Eq. (B8) can be expanded and integrated from t = to t = T. For T 
sufficiently large, all the terms vanish except the one proportional to Ak(^)^-k(^) which gives the first term on the 
right hand side of Eq. (B9). The second term on the right hand side of Eq. (|B8|) can be expanded and integrated from 



t = to t = T too. The terms proportional to Ak^k ^t^K^t^K P^^^ imaginary and give no contribution 
once added to their complex conjugate and opposite. The terms proportional to A^idtA^ and dtA^i A^, once added 
to their co mple x conjugate, both vary as S^IAkP and are easily integrated to give the last term on the right hand 
side of Eq. (|B9]). 



( dt JdxJ dyE(r,t).S^D(r,t) = 

eo5' VuK(^).u_K(^)icj^e^-(c^) / dt[Aj^{t)A_j^{t)exp{-2iujt)] 4 

K 



2.,2dK-M],^^|,^ (B9) 



Let us now calculate the energy associated with the term H.^ 



— 2]!^^ in Eq. (|B1J). 



(BIO) 



H(r, t) = —V X A(r, t) = i "S^ exp(iK.r)bK(^;)^K(*) exp{—iojt) + c.c. 
MO ^ 

dtB{r,t) = 5iVx A(r,t) = i^ exp(zK.r)bK(^) [dtAj^{t) — iujAj^{t)] ex.p{—iujt) + c.c. (Bll) 



K 



Hence 



H(r,t).StB(r,t) = V exp[i{K ^K').r]hj^{z).hj^.{z)Aj^{t) [dtAj^^{t) - iu' Aj^.{t)] 

^0 K,K' 

^ exp[z(K - K').T]hj^{z).hl^^{z)Aj^{t) [dtA^{t) + iuA^{t)] exp[-i(a; - u')t] + c.c. (B12) 



Mo 



K,K' 



Integration over the surface S as previously for E.S^D yields: 

[dx [dy H(r,t).S,B(r,t) = -^y2hj^{z).h_j^{z)Aj^{t) [S,A_K(t) - icjA_K(t)] exp(-2ia;t) 



Mo 



^|bK(^)|'AK(t) + iojA'j^it)] + C.C. (B13) 



K 



The first term on the right hand side of Eq. (B13) can be expanded and integrated from t = to t = T. For T 
sufiicie ntly l arge, the term proportional to Aj<i{t)dtA_j^{t) vani shes, and we get the first term on the right hand side 
of Eq. (B14). The second term on the right hand side of Eq. (B13) can be expanded and integrated from t = to 



t = T too. The term proportional to Ai^{t)Al^{t) is a pure imaginary, and gives no contribution once added to its 
complex conjugate. The term varying as A-K{t)dtA^(t), once added to its complex conjugate, varies as dt\Aic{t)\'^ 
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and is easily integrated to give the last term on the right hand side of Eq. (B14). 



j dt JdxJdyli{r,t).dtB{r,t) = 

VzbK(^).b_K(^) / dt [AK(t)cjA_K(t)exp(-2za;t)] 
^ Jo 



Mo 



c.c. 



f ^IbK^lVKp. (B14) 



K 



We now add Eqs. (B9) and (B14). The first terms on the right hand sides of these two equations cancel each other 
(so do their complex conjugates), so that we get 



J dt Jdxjdy [E(r,t)AD(r,t)+H(r,t)AB(r,t)] = ^eoSuj' 



duj 



UKW|^ + -^|bK(^)|^ 



Before integrating this expression over let us remark that it can be cast under a somewhat clearer form 



J dt jdxjdy [E(r,t).S,D(r,t) + H(r,t).a,B(r,t)] = 5^ 



eo 



dto 



Ek(^) 



Mo 



IBkWI^ 



(B15) 
. (B16) 



where Ek(^) = iojAj^uj^{z) and Bk(^) = ^^KbK(^) so that E(r, t) = exp(zK.r)EK(^) exp(— + c.c. and 
B(r, t) = exp(zK.r)BK(^) exp(— zcjt) + c.c. Eq. (B16) gives the energy per unit length along the z direction 



of surface plasmons. cqS 



d[u;ej(u;)] 
do; 



Ek(^)P and -^5'|Bk(^)P are then the electric and magnetic contribution of each 



mode K to the former energy per unit length. 

We now would like to integrate Eq. (B15) over z. In both half spaces z > and z < (j = 1, 2), we have |uk(^)P 

^exp[-2Im(7,)^] (l + bK(^) = {KK^jjz) x exp(i7,z)(K - ^z) = exp(z7,z)f^ x K (we 



used the property z.K = 0), hence |bK(^)P = exp[— 2Im(7j)z] 



7j 



-. Now writing 



/oo 
dz |uk(^) 
-oo 



du 



/OO 
dz \hK{z) 
-OO 



— y — 



2 4 



(B17a) 
(B17b) 



Integrating Eq. (B15) over z using Eqs. (B17), we get 
U = J^ dt jd^r [E{r,t).dt'D{r,t) + H{r,t).dtB{r,t)] = 



du 



K — J 

We now use the degree of freedom to set L{uj) as we wish, to simplify this equation. We set 



7i 



2 2 



I^Kp. (B18) 



duj 



2 2 



(B19) 



Using Eq. Q, the definition of 7j, and e{uj) < —1 at the frequencies of surface plasmon of a single interface, Eq. (B19) 
writes 



L{lo) 



2|7i| ' 4|72| L -eH ^ ^, 



-e{u:) 



(B20) 
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This is equivalent to set 



J-oo ^ 



|uk(^)| +^\hK{z)\ 



(B21) 



(the term inside the brackets comes from Eq. ( B15[ )). Using the expressions of |uk(^)P and |bK(^)P given above, 
Eq. Q and the definition of 7^, this gives a normahzation condition on uk(^) 



00 1 
dz - 



d[wej(w)] |l + e(w)| 



Eqs. {B21) and (B22) can also be written 

"eo d[uj€j{uj)] 



dz 



'1 + kHI 



1 



|uk(^)P = 1. 



2 d. .Ek(.)P + -|Bk(.)P 



with the notations used in Eq. (B16). 



With this choice of L{uj) or this normahzation condition of uk(^), Eq. (B18) can be written: 

U = Y,^oSoj'2\Aj^\'. 



(B22) 



(B23) 



(B24) 



K 



Appendix C: Quantum calculation of the surface 
plasmon emission rate of a dipole 

We follow the same steps as in the derivation of the 
photon emission rate by a two-level system in a vacuum. 
The fundamental and excited states of the two-level sys- 
tem are denoted |1) and |2) respectively, associated with 
energies Ei and E2. We define a circular frequency cjq, 
so that E2 — El = hujQ. The two- level system is initially 
in its excited state, and there are surface plasmons so 
that the global initial state can be written: \i) = |2, nx). 
In the final state, the dipole is in its fundamental state 
and a surface plasmon has been created in a mode K. 
We denote the final global state |/) = |l,nK + !)• The 
interaction Hamiltonian is — D.E, where D is the electric- 
dipole moment operator and E the quantum electric field 
operator introduced in section [IV| at the position of the 
emitter. The emission rate (inverse of the lifetime r of 
the excited state) is given by Fermi's golden rule: 



7 



27T 



^|(/|D.E|z)p(5(^2-^i-M. (CI) 

/ 



This expression can be rewritten as a sum over the modes 
K: 

7= ^^Mk(^(^2-^i-M, (C2) 



K 



where Mr = 1(1, riR + 1|D.E|2, nK)p. We note 
(2|D|l)=Di2. 



(16), we obtain the following matrix element: 



Mk = 



2eoS 



|Di2.ui,K(^)r (hk + I). (C3) 



In this equation, the tik term stands for the stimulated 
emission and the constant term 1 accounts for the spon- 
taneous emission. This section is devoted to the spon- 
taneous emission so that we do not consider the term 
associated to njo 



Using Eqs. (C2) and (C3), substituting a continuous 
sum over the vectors K in polar coordinates to the dis- 
crete sum XIk' writing jspont instead of 7 we 
get 



Ispont = ^ I dK K 



Tiuj 



(27r)2 2eQS 



5 {E2 - El- huj) 



I 

Jo 



d^ |Di2.ui,KWr. (C4) 



The integration over the directions ^ of K is performed 
usme; Eq.(|5|: 



/ d^ |Di2.Ui,k(^)|' 

Jo 



27r 



Leff{z, di2,OJo) 



(C5) 



Leff{z^di2^oJo) is defined by Eq. ( |19| ). The spontaneous 
emission rate of surface plasmon can then be cast in the 
form given by Eq. (18). 

Averaging jspont over all the possible directions (6>, (j)) 
and all the precession angles i/j of D12, we get the "total" 
spontaneous emission rate of surface plasmons 



7spont, totaz( I D12 I, CJo, z) 



Using the former expression and Eqs. (13), (14) and ^^'^ 





p27r 




[ de s'mO 


/ 


'/ 


Jo 


Jo 


Jo 



27r 



(C6) 



where R6',(/),^(Di2) is D12 rotated by Euler's angles for 
nutation, precession and intrinsic rotation ^, cj) and tjj 
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respectively. The result of the integration over 6>, (j) and where f (K, is given by 
il) is given by Eq. (23). 



Appendix D: Derivation of the emission rate of a 
dipole in the classical lossy case 

The aim of this section is to derive an explicit form of 
the Purcell factor due to the presence of surface plasmons 
by using the Green's tensor approach. Using Ref. 35 , one 
can write the surface plasmon contribution to the Green's 
tensor evaluated at the position of the source r: 



,(r,r,cj) = j 



d^K 

(2^)2 



f(K,z,z) f*(-K,z,z) 



UJ ^ UJ 



sp 



^ ^ K ^ K 

f (K,2;,z) = -g{K,uJsp){^ z)(K z) exp(2z7i2;) 

7i 7i 

(D2) 

^^4^^^(^>w«p) and 



as z > 0, with g{K,ujsp) 



y^^) Injecting Eq. in Eq. (24), we obtain: 



^P,d(di2,cJ,z) = 



Gttc 



/ 



(27r)2 



Im 



dj2. f (K, Z, Z)di2 dl^.f * (-K, Z, Z)di2 



•^sp 



ijj -\- (jj 



sp 



(D3) 



Writing FK(di2,w,z) = dJa- f (K, ^)di2 = Fi^{du,oj, z) + iF^{du,oj, z) and p'^^^{uj) = (^_^, j2+^„2 , 
7 — 772 , we find: 



Fp,ci{di2,uJ,z) 



6ttc 



d^K 

(27r)2 



-F^(di2,w,2y:„. H + F^(di2,w,^)p'_„. H 



• (D4) 



This expression is to be compared to Eq. (C2) in the non lossy case. Both expressions are written as sums over the 
modes K, of the contribution of ea ch m ode to the spontaneous emission rate. When losses are low, the last two 
terms on the right hand side of Eq. (D4) are antiresonant, as p'^^* (cj) and p'_^* {u) are centered around —uj^p. The 
second term has also a small contribution as the average value of p'^^^^iuj) is 0. The main contribution comes from the 
first term, as p'^ {u) goes to —iTi8{uj — ujsp) in the non lossy limit. Its extremum value is ^77^. Starting again from 



Eq. (Dl ), we write 



G5p(r,r,cj) 



1 



(2^)^ 70 



/ dKK \ / d^ f(K,z,z) -^-^ f*(-K,z,z) 



aiming at performing the integral over the directions ^ of K. We first calculate: 



/.27r 




/ dO 


'f(K,z,z)" 


'0 



where \J{z^uJsp) = ex.p{2ijiz) [^(xx + yy) — e{uJsp)zz]. Inserting Eq. ( |D6[ ) in Eq. ( |D5[ ) yields: 



^ 1 r 

G,p(r,r,a;) = ^ 



dKK 



- U [z, uJsp ) H ^ , , ^ — ^ (^5 ^sp ) 



Inserting Eq. ( |D7| ) in Eq. ( |24[ ), one finds: 

3c 

^p,d(di2,cJo,^) = — / di^K Im 
^0 Jo 



^0 ~ ^sp 



di2.U(z,a;5p)di2 



UJ -\- UJ 



di2.U*(z,a;5p)di2 



(D5) 



(D6) 



(D7) 



(D8) 
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with d*2.U(2;, Wsp)di2 = exp{2i^iz) [^\<ii2,//\'^ - e{i^sp)\di2,z\^]- In the non lossy case, ^(i^, w^p) and dj2-U(z, Wgp)di2 
are real, and Im-^j^ — goes to'*-'- —inS{uj — LOgp)- Hence, assuming w > and using the expressions of dj2-U(i;, Wsp)di2 
and g{K,uJsp) given above, one gets Eq. (25). We need the relation 



1 



(D9) 



in order to prove that Fp(di2,cJo, z) = Fp,cZ (di2, cjq, z). Using the definitions of R{K^ujsp) and of given above, we 
get 



71 



de 
dcj 



de 
da; 



272 



Inserting ^ 



^ ~ ^] ( jDlQh and using 71 = i|7i|, 72 = -i|72|, we find 



l7i 



1 w/c^A dwe 



w 2I72I/ dw 



Using Eq. Q and the definition of 7^, we can simplify the terms inside the parentheses and write 



euj/c^ e eui/c^ w/c^l — edwe 
|7i| |l + e| ~ 2I72I " 2^rT^ d^7 



Multiplying this expression by 



we get 



e+ 1 



1 1 — e dwe 



|7i| 2I72I 2I72I e dw 



Comparing this expression to Eq. ( |B20[ ), we can write 

jR-\K,LUsp) = 2L{lu) 



U! 



(DIO) 



(Dll) 



(D12) 



(D13) 



(D14) 



which gives Eq. (D9). 
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